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ABSTRACT 

An explicit solution of the linearized equations of motion is used 
to obtain the deviations in end conditions, Ar and A0, and the addi- 
tional burning time, At, expressed as integrals of functions of thrust 
angle deviations, AX. Under the constraint that Ar = A0 = 0, the cal- 
culus of variations is applied to minimize At. The resulting Euler- 
Lagrange equation evaluated at AX = 0 gives a simple relationship which 
is used to show that At is a second order function of AX. This function 
is evaluated numerically for the second stage of an early SA-6 design. 
Results show that, if the mission is accomplished, the thrust angle may 
differ by as much as 2 degrees throughout the second stage with a pro- 
pellant penalty of only 34 pounds. It follows that the capability of 
meeting the mission is the prime requisite of a guidance function and 
optimality is then only a second order consideration. 
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DEFINITION OF SYMBOLS 


Symbol 

t 

t Q 

W 

Ar 

At 

At 0 

ZM 

A0 

AX 

X 


Definition 

time measured on a standard trajectory 

second stage ignition time on the standard trajectory 

second stage cutoff time on the standard trajectory 

propellant flow rate 

radius error at cutoff 

deviation in second stage burning time 

deviation in second stage ignition time from standard 

deviation in propellant consumption from standard 

error in angle between the velocity vector and the posi- 
tion vector at cutoff 

X(t + At 0 ) - Xg(t) 

thrust angle 

thrust angle on the standard trajectory 
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TECHNICAL MEMORANDUM X-53249 - 


NEAR-OPTIMUM GUIDANCE - AN ANALYSIS OF FUEL PENALTY 


SUMMARY 

An explicit solution of the linearized equations of motion is used 
to obtain the deviations in end conditions, Ar and A9 , and the addi- 
tional burning time. At, expressed as integrals of functions of thrust 
angle deviations, AX, Under the constraint thatAr = A0 = 0, the cal- 
culus of variations is applied to minimize At. The resulting Euler- 
Lagrange equation evaluated at AX = 0 gives a simple relationship which 
is used to show that At is a second order function of AX. This function 
is evaluated numerically for the second stage of an early SA-6 design. 
Results show that, if the mission is accomplished, the thrust angle may 
differ by as much as 2 degrees throughout the second stage with a pro- 
pellant penalty of only 34 pounds. It follows that the capability of 
meeting the mission is the prime requisite of a guidance function and 
optimality is then only a second order consideration. 


I. INTRODUCTION 

The task of constructing a guidance function to determine the thrust 
angle, X, that meets a given mission under the constraint of minimum pro- 
pellant consumption can be considerably simplified by first investigating 
the penalty for non-optimality. It is assumed that minimum propellant 
consumption is equivalent to minimum burning time, and it will be shown 
that any additional burning time required is a second order function of 
thrust angle deviations from optimum provided the mission is met. A 
numerical example is included which covers the second stage of a 100 n.m. 
orbital mission for an early SA-6 design. The results show that any 
guidance function which meets the required end conditions is acceptable 
even though it differs from the calculus of variations solution by as 
much as one or two degrees throughout the entire second stage. 


II. EULER- LAGRANGE EQUATIONS 


It is assumed that the calculus of variations solution is known 
for a given trajectory. If the thrust angle, X, deviates slightly 
from optimum, it has been shown in Reference 1 that the errors in end 
conditions Ar and A0 and the additional burning time At can be adequately 
expressed as follows: 


Ar 


f(AX,t) dt. 


(1) 


t 

r 

j 


t 

o 


n 

A0 = J' g(AX,t) dt, 
t 


( 2 ) 


n 

At = J" h(AX,t) dt, 
t 


(3) 


where 


f (AX, t) = f ! ( t ) AX. + f 2 (t) AX 2 + ... , (4) 

g (AX, t) = gl (t) AX + g 2 (t) AX 2 + . . . , (5) 

h(AX,t) = hjCt) AX + h 2 (t) AX 2 + ... . (6) 

Minimizing equation (3) under the constraint that equations (1) 
and (2) satisfy prescribed conditions leads to the Euler-Lagrange equa- 
tion [2, p. 51] 


Bf* _d_ BjT 

Bax " dt ba£ 

where 


f" = h + A x f + A 2 g. 


(7) 


( 8 ) 


2 



Since f' v is independent of AX, equation (7) degenerates to the following: 


Sax Sax Al Sax 


A 



which yields the following necessary condition; 


(ha. + Mi 


+ A 2 g x ) + 2 (h 2 + Axf 2 + Aog 2 ) AX + ... - 0, 


(9) 


Since the standard was assumed to be a calculus of variations solution, 
equation (9) is satisfied for AX = 0 and the following relationship is 
determined: 


h-i "b A if i "b A 2 g i — 0 • 


( 10 ) 


Table I shows the numerical values obtained for h 1? f ls 
following values of A a. and A 2 were obtained by the method of 


and g ! . The 
least squares. 


Aa. - -.08734 sec/km 
A 2 = .8589 sec/ deg. 


These values satisfy equation (10) within the numerical accuracy with 
which h l5 f l5 and g ± were determined. This is illustrated in Table I 
where the quantity h x + A x fi + A 2 g! is shown. 


III. SECOND ORDER EFFECT 
Fr om e qua t ion (10). it follows tha t 
h x AX + Aifa/NX + A 2 g! AX = 0 
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TABLE I 


(10“ 2 km/deg sec) 
-3.680 
-3.738 
-3.772 
- 3.774 
-3.738 
-3.654 
-3.507 
-3.280 
-2.944 
-2.458 
-1.758 
- .727 

= -.08734 sec/km 
A 2 - .8589 sec/deg 


Si -4 

(10 2 /sec) (10 2 /deg) 
.0454 -.3610 

.0536 -.3730 

.0628 -.3836 

.0732 -.3927 

.0850 -.3996 

.0986 -.4038 

.1144 -.4044 

.1330 -.4005 

.1554 -.3903 

.1828 -.3715 

.2177 -.3405 

.2643 -.2907 


h 2 + Anf i + A 2 gi 
(10 -2 / deg) 

-.0006 

-.0005 

-.0002 

-.0002 

-.0001 

.0000 

.0002 

.0002 

.0003 

.0002 

.0000 

-.0002 


and 


t 

n 


t 

n 


t 

n 



x AX dt = -A x 


f x AX d t — A 2 


g 1 AX. d t , 


t 

o 


t 

o 


( 11 ) 


From equations (4) and (5) together with equations (1) and (2), 
the following relationships are determined: 


t 

n 


t 

n 


- A x 


f. x £X dt = -A x Ar + A x 


[f^AX 2 + ... ] d t , 


t 

o 


t 

o 


and 


A 2 


t 

n 


g x AX dt = 


t 

o 


t 

n 

r 

A 2 A0 + A 2 / [ga^X 2 + ...] dt. 
u 
t 

o 


Substituting these values into equation (11) yields 


t 

n 


h-j_AX dt = - AjAr - A^G + 


t 

n 


[ (A x f 2 + A 2 g 2 ) AX 2 + 


] dt. 


t 

o 


t 

o 


This, together with equations (3) and (6), gives the following expression 
for At. 


At = -A x Ar - A^AQ + / [ (h 2 + A x f 2 + A 2 g 2 ) AX + . . . ] dt. 


( 12 ) 
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If the required end conditions are met, Ar = AG = 0. Then, 


At 


t 

n 


[(h 2 + A x f 2 + A 2 g 2 ) + 


] dt. 


t 

o 


( 13 ) 


and the additional burning time At resulting from non-optimum X is a 
second order function of AX. The sufficient condition that t is a 
local minimum is that h 2 + A x f 2 + A 2 g 2 > 0, t Q ^ t ^ t n . 


IV, NUMERICAL RESULTS 


Table II shows the numerical values obtained for h 2 , f 2 , g 2 arid 
(h 2 + A 2 fe + A 2 g 2 ). In addition, the quantity h*(t-[) = (h 2 -f A 2 f 2 + A 2 g 2 ) At^ 
is shown in Table III. This quantity can be used to obtain a good approxi- 
mation to the integral of the second ordered term in equation (13) by the 
following summation. 

12 

At = y h*(t.) AX 2 ( t i ) . (14) 

i=l 


This expression, in turn, has the following upper bound: 

12 ^ 

At S AX 2 ) h*(t.), (15) 

i=l 

which, for this example, is 

At A .0409 AX 2 , (16) 

m 

where A^ is the maximum value of |AX(t)|, t Q S t A t n . 

The propellant loss AW can be expressed as follows: 

Atf = WAt. 

For this example, W = 208 lbs/sec, which yields the following bound: 

AW A 8.51 AX 2 (lbs/deg 2 ). 
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(17) 



TABLE II 


t f 2 g 2 

(sec) (10“ 3 km/deg 2 sec) (10 -o /deg sec) 


h 2 (h 2 + Aif 2 + A 2 g 2 ) 

(10" 4 /deg 2 ) (10" 4 / deg 2 ) 


160 

-.3831 

.8899 

.2168 

.6278 

200 

-.3516 

.8798 

.2606 

.6433 

240 

-.3194 

.8663 

.3094 

.6628 

280 

-.2863 

.8482 

.3641 

.6870 

320 

-.2524 

.8239 

.4260 

.7172 

360 

-.2175 

.7915 

.4969 

.7548 

400 

-.1819 

.7483 

.5789 

.8021 

440 

-.1458 

.6907 

.6755 

.8622 

480 

-.1096 

.6131 

.7915 

.9399 

520 

-.0742 

.5074 

.9344 

1.0428 

560 

-.0408 

.3602 

1.1167 

1.1832 

600 

-.0118 

.1485 

1.3604 

1.3834 


A x = -.08734 

A 2 = .8589 sec/deg 


4 


7 




V. CONCLUSIONS 


Equation (17) shows that the propellant penalty from non-optimum 
guidance is less than 34 pounds even if the thrust angle deviates from 
optimum by qs much as 2 degrees throughout the second stage. This 
clearly indicates that, if a guidance function can be determined which 
meets the mission within satisfactory limits and predicts X within one 
or two degrees of the optimum value, there is very little to be gained 
by expending much effort to further reduce the prediction error in X. 
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